Thermal rho and sigma mesons from chiral symmetry and unitarity 
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We study the temperature evolution of the p and a mass and width, using a unitary chiral ap- 
proach. The one-loop nn scattering amplitude in Chiral Perturbation Theory at T 7^ is unitarized 
via the Inverse Amplitude Method. Our results predict a clear increase with T of both the p and a 
widths. The masses decrease slightly for high T, while the pnn coupling increases. The p behavior 
seems to be favored by experimental results. In the a case, it signals chiral symmetry restoration. 
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One of the outstanding phenomena related to heavy 
ion collisions is the flatness of the dilepton spectrum near 
the mass of the p meson, which is so clearly visible in 
many processes involving hadrons and electromagnetic 
probes. This flatness has been observed by the HELIOS 
and CERES collaborations J|, |j and has been the sub- 
ject of widespread discussion. Dileptons and photons pro- 
vide neat signals of the early stages of the quark-gluon 
plasma and its subsequent evolution into a hadron gas 
H . In fact, the most credible explanation of the absence 
of a prominent hill in the dilepton spectrum is a change 
in the mass and width of the p due to its interactions 
with the hot hadron gas 0, EL 0, Since the baryons, 
with a large forward momentum, have almost escaped 
the central collision region, this gas is composed mainly 
of pions. Our aim is to study the thermal evolution of 
the p mass M p and width T pi from the first principles of 
chiral symmetry and unitarity in irir scattering. 

What happens to the p in extreme conditions is a 
hadronic physics problem, involving non perturbative 
physics and hence difficult to be treated. Prior to this 
work, a copious number of models and estimations have 
appeared. In most of them T p increases with tempera- 
ture, simply as a consequence of stimulated emission in 
the pion thermal bath or, equivalently, because the effec- 
tive phase space increases fi §. This behavior is often 
interpreted as a deconfming effect, or hadron "melting". 
As for the mass, Vector Meson Dominance (VMD) im- 
plies that M p changes very little at low temperatures 
|^|, 0. As T approaches the critical temperature, 
earlier works claimed that M p increases [||, |ll|, [l^] but 
the analysis of experimental dilepton data seems to favor 
a decreasing behavior ^ [7J. Let us remark that in all 
these works, the p is introduced as an explicit degree of 
freedom and often a dilute pion gas is assumed, so that 
the thermal effects appear, to leading order, only through 
the pion distribution function and not through the inter- 
action details. Other approaches include the NJL model 
plU , where M p and T p slightly decrease (but there is an 
spurious quark threshold near M p ) as well as qq wave- 
function analysis in the tt channel yielding a decreasing 

width in. 



In this work we will use a thermal treatment of the 
effective degrees of freedom, the pions in the aftermath 
of the collision at moderate temperatures. The guiding 
fundamental principles will be just chiral symmetry and 
unitarity. We will build on a previous work [flEf where 
the T 7^ 7T7T scattering amplitude has been calculated 
to one loop in Chiral Perturbation Theory (ChPT). De- 
manding unitarity, we will construct a non-perturbative 
amplitude reproducing the expected behavior for ther- 
mal resonances. Our amplitude has the correct analytic 
structure, without spurious cuts, and resonances are not 
introduced by hand. 

The most general framework comprising the QCD chi- 
ral symmetry breaking pattern is ChPT Jr], [ItJ , see (ll| 
for reviews, where observables are calculated as expan- 
sions in p/(47r/ ff ), p denoting any pion energy scale (in- 
cluding the temperature) and f„ ~ 92.4 MeV. Despite 
its success, ChPT is limited to low energies (usually, less 
than 500 MeV) and low temperatures and it is not able 
to generate resonances. Thus, over the last few years, 
there has been a growing interest to extend the ChPT 
applicability range to higher energies and to reproduce 
resonances within a unitary chiral approach, which we 
briefly review. At T = 0, unitarity for the S'-matrix 
(S'S = 1) implies the following relation for partial waves 



lma u (s) = a(s)\au(s)\ r/: 



(1) 



for s > 4m 2 and below other inelastic thresholds, where 
<r(s) = yl — 4m*Js is the two-pion phase space and au 
denotes the projection of the 7T7t elastic amplitude with 
isospin / and total angular momentum J in the center of 
mass frame. Eq.(|l]) is only satisfied perturbatively within 
ChPT, i.e, if we write the perturbative series for any par- 
tial wave as a = 0,2 + a a + • • ■ where is 0(p k ), then 
one has Im 02 = 0, Im = aa\ and so on. Hence, devia- 
tions from eq.(Q) are more severe at high energies, and in 
particular near the resonance region, where the bounds 
imposed by unitarity are saturated. The ChPT series, 
which essentially behaves as a polynomial, is unbounded 
and cannot reproduce resonances, which show up as poles 
of the amplitude in the complex plane. 
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In fact, from eq.([j]), any partial wave should satisfy 
a = l/(Rea _1 — ia) on the real axis below inelastic 
thresholds. A unitarization method is just one way of 
approximating Re a -1 , thus introducing some model de- 
pendency, but since we want to ensure chiral symmetry, 
and the correct low energy behavior at T = 0, we will 
use the one-loop ChPT result. This is called the IAM 
at T = 0, which can be recast as a IAM — 0,2/(0,2 — a±) 
p9[ . The single channel IAM amplitude satisfies eq.(|l|) 
exactly and at low energies it follows the ChPT result up 
to one loop. In addition, the IAM reproduces the scatter- 
ing data for real energies above the two pion threshold up 
to 1 GeV, where the elastic approximation breaks down, 
and it can be continued into the complex s-plane, yield- 
ing correct a and p poles in the second Riemann sheet. 
We point out that the IAM is nothing but the [1,1] Padc 
approximant of the ChPT series in squared energy, mass, 
or temperature over 

As long as they contain the 0(p 4 ) tree level terms, 
other chiral unitary approximations, both for SU(2) or 
SU(3) ChPT, either based on the IAM with coupled 
channels p0| , or the IAM with higher orders pi]], or in- 
spired in Lippmann-Schwinger or Bethe-Salpeter equa- 
tions p2| , or mixed formalisms |23| yield equivalent re- 
sults for the p and a channels. In particular, they repro- 
duce the experimental phase shifts with compatible sets 
of chiral parameters, and they generate poles associated 
to the a and p resonances whose position in the second 
Riemann sheet agrees for all the above mentioned meth- 
ods, which therefore describe resonances with the same 
masses and widths. These unitarized approaches also al- 
low to study finite baryon density effects on the change 
of the sigma properties in the nuclear medium 24 that 
suggested a decrease on both the sigma mass and width 
as the nuclear density increases. As a consequence of 
these effects, it is expected j25) a shift of strength of the 
two pion invariant mass distribution in 7 A — > A7r°7r°, 
which has been recently confirmed experimentally ]p6| . 
In particular, using a chiral unitary approach, this shift 
is interpreted as an in medium modification of the a pole 
towards lower masses and widths Nevertheless, since in 
this work we are interested in the p and a mesons it is 
enough to work with the sing le channel IAM to 0(p i ) 
that we have just described. 

Back to T 7^ 0, the thermal amplitude can be defined 
by considering T = initial and final asymptotic states 
and calculating the T ^ four-pion Green's function 
Jl5| . To one loop in ChPT, and in the nn c.o.m. frame 
(at rest with the thermal bath) it satisfies the perturbative 
unitarity relation jlBf : 



where 



lma 4 (s;T) = a T (s) [a 2 (s)Y 



a T {s) = cr(s) [1 + 2n B (v^/2) 



(2) 



(3) 



is the thermal phase space and Ub{x) = (exp(x/T) — 1) _1 
is the Bose-Einstein distribution function. Recall that 
the lowest order a 2 is T-independent. 

Therefore, the natural unitarized version of the ther- 
mal amplitude in ChPT should be: 



JAM 



a|(s) 



a 2 (s) - a 4 (s;T) 
which satisfies the exact elastic unitarity condition 



Imo 



IAM 



( S ;T) = a T (s)\a IAM (s;T) 



(4) 



(5) 



and reproduces the low energy results of (thermal) ChPT 
in p5| . Besides, as we will see below, it has the proper 
analytical behavior and, for the appropriate values of the 
chiral parameters, it is able to reproduce resonances like 
the p as poles in the second Riemann sheet. 

Some remarks are in order here: We are assuming that 
the exact thermal version of eq.([l]) holds, feature repro- 
duced by the IAM in eq.(||). This assumption will prove 
to be reasonable in view of the results shown below. Nev- 
ertheless, it is important to remark that such assump- 
tion implies in particular that only two-pion states are 
available in the thermal bath. This is equivalent to a 
dilute gas approximation. In other words, the Ub term 
in eq.@ must remain small compared to one so that we 
can neglect higher orders in density like 0{n 2 B ) which 
would spoil the simple algebraic unitarity relation given 
by eq. (||) |l5| . This implies, for instance that p — tt scat- 
tering, which in our approach is regarded as a three-pion 
effect, would be suppressed by the low density. Note 
that, alternatively, we can view eq. (Q) also as the [1,1] 
Pade approximant of ChPT when counting the powers of 
momenta, masses or temperature over 1//^, since T is 
0(p) in the chiral expansion. Again, this counting would 
be spoiled for large ns(yi), which typically weights the 
thermal corrections. Let us finally remark that the IAM 
has been extended to deal with other intermediate states, 
describing successfully the T = data in all meson-meson 
channels up to 1.2 GeV, within a coupled channel formal- 
ism pof . In such case, the amplitudes satisfy a matrix 
version of the unitarity relation in eq.(|l]). This would 
be the natural extension of our approach in order to deal 
with other intermediate coupled states, like K or rj which 
could be relevant at high temperatures f27|| . However, 
it would require a thermal generalization of the matrix 
unitarity relation and the one-loop calculation of the ad- 
ditional coupled amplitudes, which lie beyond the scope 
of this work. 

Before proceeding to the detailed calculation of the 
IAM thermal amplitude in eq.(^), we will provide a sim- 
ple argument as to why our method can actually give rise 
to the expected thermal behavior for the p width. As it 
is well known, in most cases (like the p) a resonant be- 
havior can be reproduced on the real axis by means of a 
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Breit-Wigner parameterization of the partial waves: 



BW 



(s;T) 



Rt{s) 



M 2 



iT T M 7 



(6) 



where Mt and Yt are the thermal mass and width of the 
resonance Q and Rt(s) is a smooth real function near 
s = M T , which can be related to the pmr vertex (see 
below). The parameterization in eq.(||) applies only for 
s ~ M T and for narrow resonances (IV <C Mr). Com- 
paring eq.(Q) with eq.@ at s = M T one readily gets 
Rea4(M T ) = <i2(M T ) (the resonance mass condition) 
and, using eq.(|), T T M T = -R T (M T )a T (M T ). There- 
fore, assuming that the thermal corrections to Rt and 
to Mt are much smaller than those to Tt, i-e, Rt — Ra 
and Mt — Mq we would get 



r [1 + 2n B (M /2)] 



(7) 



Hence, in this limit the thermal IAM yields an increas- 
ing resonance width driven only by the available thermal 
phase space eq. @ for a p at rest || || . The above result 
takes into account the stimulated emission p-^ kit and 
absorption tttt — ► p from the thermal bath |L5| and gives 
the dominant effect at very low temperatures, as our full 
analysis below confirms. This approximation indicates 
that the unitarity requirements on the amplitude capture 
the qualitative thermal resonance behavior. Note that, 
from the resonance mass condition, taking Mr — Mq is 
equivalent to ignoring the T-dependence in Rectus;! 1 ). 

Therefore, by using the full thermal amplitude <24(s; T) 
in |L5|, we will calculate below both the Mr corrections 
and the deviations from eq. (Q) . Moreover we will find the 
analytic continuation of the amplitude to the complex 
plane, so that we can describe the resonances as poles 
of the thermal amplitude. This is particularly important 
for the a, whose description in terms of eq.(^) is not so 
appropriate due to its large width. 

Let us note that the breaking of Lorentz invariance 
of the thermal formalism allows for a definition of a 
"transversal" and a "longitudinal" mass, however, in our 
case, since we are working in the c.o.m. frame, where the 
p is at rest with the thermal bath, both mass definitions 
coincide ]jT| . 

In the c.o.m. frame, the thermal one-loop amplitude 
can be written in terms of the loop functions |32] : 



AJJ( S ;T) =-- 



AE 2 



AJ* u (i;T) 
A J? (i;T) 



1 



dq 



gn B {E q ) 
E„ 



log 



2q + 4~t 



dqqE q n B {E q ) log 



v /= i-2g 
2q + V^i 



V=i-2q 



for real s > Ami an d real t < 0, where AF(T) = F(T) 



F(0), E\ 2 = q 2 + ml, J^(s;0) is given in E7) (after the 
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FIG. 1: I — J — 1 phase shift for different temperatures. For 
the data see and references therein. 
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FIG. 2: Temperature evolution of the p mass, width and pivn 
coupling. The dashed line corresponds to eq.(Q), the dotted 
line to the real axis IAM, eq.(^), and the solid line is the IAM 
pole position. The Mq, To, go values are given in the text. 



standard MS — 1 renormalization) and Jo^^O) can be 
written in terms of JJ(t;0). Note that on the real axis 
the only imaginary part comes from Im J§(s + ie;T) = 
ctt(s)/167t 2 for s > 4m 2 , thus ensuring eq.(||) [f33[ . 

We have calculated the IAM amplitude, eq.(|J), from 
the one-loop thermal 04 (s; T). The phase shifts for differ- 
ent temperatures are shown for the p channel I = J = 1 
in Fig.jl. Note the excellent agreement with scattering 
data at T = 0, where we have fitted the SU (2) low-energy 
constants in the aoo, an, 020 channels, yielding the fol- 
lowing values for the standard dimensionless and scale 
independent SU(2) chiral parameters defined in [Jl?], 
li = —0.3, I2 = 5.6, 73 = 3.4 and I4 = 4.3, compatible 
with the recent determination in |p8| . The p mass and 
width can now be estimated using Sn(M p ) = 90° and 
T p (s) ~ M p (l - s/M 2 )tan<5n near s = M 2 @. Thus 
we obtain M = 770 MeV and T = 159 MeV. All fi- 
nite T results are now predictions. As T increases, Ft 
grows, as shown in Fig.^. The curves are shown only be- 
low the validity limit of our approach which naively is set 
by 2n B (M p /2) < 1 yielding T < 300 MeV. We remark 
that the validity of one-loop SU (2) ChPT has been es- 
timated to reach about T ~ 150 MeV but with the 
unitarization methods we are able to reach higher tem- 
peratures, as long as the density factors remain small 
(see our previous discussion). To lie on the conserva- 
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tive side, we are only showing results up to 200 MeV, 
where 2riB(M p /2) ~ 0.3. Let us still note that devia- 
tions from the naive phase space correction in eq.(^) are 
clearly sizable already at T ~ 100 MeV, precisely when 
thermal effects start being significant, the full calcula- 
tion giving a higher value for the width than eq.(Q). The 
mass changes little up to T ~ 200 MeV, consistently 
with previous analysis |7], ^, [l0[ [l^, [D| O. ft grows 
slightly up to T ~ 100 MeV {M wo ~ 775.5 MeV) and 
then decreases for higher T. In addition, in the narrow 
resonance approximation (which becomes less reliable as 
T increases) we have Rt = g\ (4m^ — M|,) /487r, gx be- 
ing the effective coupling in the VMD p7T7r vertex [|[ ||] 
with a thermal p (go ~ 6.2). Therefore, from the IAM 
Tt,Mt we find the behavior of qt plotted in Fig.|| At 
low T, gr^go (ff5o/so — 0.9991) in agreement with the 
chiral low-T analysis in || and it grows for higher T. The 



corrections are more important at finite density |29j . 

Although the direct experimental measurement is the 
dilepton spectrum, tttt scattering is still a very interesting 
process since it is strongly constrained by unitarity. In 
particular this provides relevant information about the p 
pole position, which has to be a common feature for all 
other processes where the p resonance appears. That is 
why we now turn to study the analytic continuation of 
the amplitude to the complex plane. The analytic contin- 
uation of the T = Jq is straightforward |l7], [1^] . How- 
ever, due to the loss of Lorentz covariance in the thermal 
bath, we need the analytic continuations of eqs.(||) which 
are somewhat more subtle. Since AJq(s;T) is already 
written as an analytic function for Im s 0, the same 
expression is straightforwardly continued to the complex 
plane. However, for the others we find: 



A ± J t u (t;T) 



4ir 2 V~t {Jo 



dq 



qn B (E q ) 
E a 



log 



2q 



t-2q 



±z^Tlog(l-e- fl ( t )/ T )| 



(9) 



A ± J* u (t;T) 



47r 2 ^ Uo 



dqqE q n B (E q ) log 



2q + ^/~t 
^f~t-2q 



± inT 



R 



'(i) log (l - e - R ^/ T ) 



2TR(t)U 2 (e-«W/ T ) - 2T 2 Li 3 (V^ 4 ^)] } 



(10) 



I 

II 



where R(t) — \J m\ — t/4, A + (~) denote the analytic 
continuation for Imi > 0(< 0) and Li„(z) is the poly- 
logarithmic function, analytic except for a branch cut for 
real z > 1 Q. It is not difficult to check that AJg* 2 (i; T) 
coincide with eq.(^j) at t±ie with real t < and they have 
a branch cut only for real t > 4m 2 . Thus, as it happened 
for T = 0, both 04 and a IAM have a right (unitarity) cut 
for real s > 4m 2 and a left cut for s < coming, respec- 
tively, from AJ S and AJ^ 2 g). Finally, using eq.(|), g 
the analytic continuation of the amplitude a 11 into the 
second Riemann sheet across the right cut is given by 
a n (s; T) = a IAM (s; T)/[l - 2ia T (s)a IAM (s; T)}. 

For I = J = 1, we find the pole corresponding to 
the p resonance. Its position on the complex plane as a 
function of T is shown in Fig.[| Let us recall that the 
definition of the pole position in terms of the resonance 
mass and width is s po i e — [M — iT/2) 2 , which coincides 
with the pole of eq.(|^) for a narrow resonance. In par- 
ticular, at T = 0, we have M = 755 MeV and T = 
152 MeV. The results are also plotted in Fig.[| where we 
see that the evolution of the pole mass and width agrees 
with our previous real axis calculation. For I = J = 0, 
the observed pole corresponds to the a and is plotted 
in Fig.|| too. The width also increases, essentially by 
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the increase of phase space and M a (T) decreases with 
T, as expected from chiral symmetry restoration [pl| . 
Once again, the applicability of our approach is limited 
by 2n B {M a /2) ~ 1, i.e, T < 180 MeV. 

The main conclusions of this work are the following. 
We have shown, using only chiral symmetry and unitar- 
ity, that the thermal width of the p and a mesons at 
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rest with the thermal bath grow with temperature, while 
their thermal masses decrease slightly. They can be read 
off from the real and imaginary parts of the pole posi- 
tion of the thermal tttt elastic scattering amplitude in the 
corresponding channels. For that purpose, we have uni- 
tarized and calculated the analytic continuation to the 
complex plane of the amplitude on the real axis above 
threshold analyzed in |^5|. For the case of the p we 
have also estimated its thermal mass, width and effective 
7T7T coupling from the unitarized amplitude in the real 
axis. At low temperatures, the thermal widths increase 
slightly according to the thermal phase space, while the 
masses and the effective pmr vertex remain almost con- 
stant. For higher T, our analysis gives sizable decreasing 
mass corrections, an increasing effective vertex, as well 
as significant deviations from the phase space contribu- 
tion, yielding higher thermal widths. The a mass shows 
a decreasing behavior compatible with chiral symmetry 
restoration. Our results agree with recent theoretical and 
experimental analysis, up to temperatures of 250 MeV 
and they shed light on the dilepton spectrum problem in 
Relativistic Heavy Ion Collisions. 
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